Journal: Nucleic Acids Research
Article Title: Separate transcription and splicing gene networks are linked and coordinated by the pRb–E2F pathway
doi: 10.1093/nar/gkag016
Figure Lengend Snippet: The pRb–E2F pathway regulates RNA splicing of E2F target genes. ( A ) Differential changes in splicing between WT and E2F1 Cr HCT116 cells, treated with DMSO or 1 μM T1-44 for 48 h are displayed as a heatmap of delta PSI values (ΔΨ, PSI) for all significant differential splicing events (FDR < 0.05, ΔΨ > 0.1). Each column of the heatmap represents delta PSI values of one splice event across all samples, as compared to the WT E2F1 HCT116 cells treated with DMSO (blue: reduced inclusion; red: increased inclusion). Data clustering used the Canberra distance method. Black boxes indicate splice events that uniquely occur upon T1-44 treatment, only in the presence of WT E2F1. These data were generated from three independent biological samples. Venn diagram showing the overlap between statistically significant differential splicing events (FDR < 0.05) (AS) in each treatment, as compared to WT E2F1 HCT116 cells treated with DMSO. These data were generated from three independent biological samples. ( C ) A representative immunoblot displaying input protein levels of E2F1 and symmetric dimethylation (SDMe). Actin served as a loading control. ( D ) The bar chart displays the breakdown of statistically significant (FDR < 0.05) differential splicing events observed in each of the indicated treatments, as compared to WT E2F1 HCT116 cells treated with DMSO. SE, skipped/cassette exon; RI, retained intron; MXE, mutually exclusive exons; A5SS, alternative 5′ splice site; A3SS, alternative 3′ splice. These data were derived from the analysis in Fig. and . ( E ) Annotation of genes which undergo splice events that uniquely occur upon T1-44 treatment in the presence of WT E2F1 (see Fig. ). GO biological process (GO:BP) and Reactome gene sets were used for pathway analysis in Metascape. Biological terms connected with cell cycle, stress responses, and DNA damage are highlighted in red. The number of genes enriched in each category is displayed to the right of each bar. ( F ) Differential changes in splicing between WT and Rb Cr MCF7 cells, treated with DMSO or 1 μM T1-44 for 48 h are displayed as a heatmap of delta PSI values (ΔΨ, PSI) for all significant differential splicing events (FDR < 0.05, ΔΨ > 0.1). Each column of the heatmap represents delta PSI values of one splice event across all samples, as compared to the WT Rb MCF7 cells treated with DMSO (blue: reduced inclusion; red: increased inclusion). Data clustering used the Canberra distance method. Black boxes indicate splice events that uniquely occur upon T1-44 treatment, only in the presence of WT Rb. These data were generated from three independent biological samples. Venn diagram showing the overlap between statistically significant differential splicing events (FDR < 0.05) (AS) in each treatment, as compared to WT Rb MCF7 cells treated with DMSO. These data were generated from three independent biological samples. ( H ) A representative immunoblot displaying input protein levels of Rb and SDMe. GAPDH served as a loading control. ( I ) The bar chart displays the breakdown of statistically significant (FDR < 0.05) differential splicing events observed in each of the indicated treatments, as compared to WT Rb MCF7 cells treated with DMSO. SE, skipped/cassette exon; RI, retained intron; MXE, mutually exclusive exons; A5SS, alternative 5′ splice site; A3SS, alternative 3′ splice. These data were derived from the analysis in Fig. and . ( J ) Annotation of genes which undergo splice events that uniquely occur upon T1-44 treatment in the presence of WT Rb (see Fig. ). GO biological process (GO:BP) and Reactome gene sets were used for pathway analysis in Metascape. Biological terms connected with cell cycle, stress responses, and DNA damage are highlighted in red. The number of genes enriched in each category is displayed to the right of each bar. ( K ) WT E2F1 and E2F1 Cr HCT116 cells treated for 48 h with 1 μM T1-44 or DMSO as indicated. An RT-PCR was performed to measure the inclusion of VCAN exon 7, MDM1 exon 4, METTL6 exon 3, or REV3L exon 3 in RNA transcripts from the cells. Displayed is the mean inclusion/exclusion ratio, with SD. Significance was calculated by ANOVA using Sidak’s multiple comparisons test. A diagram indicating the exon (boxed in grey) and intron (black lines) structure of each gene around the skipped exon (boxed in yellow) of interest is included. The splicing that gives rise to the exon included and excluded transcripts is also displayed, with specific primer pairs used in QPCR shown as blue arrows. A representative immunoblot is included to display input protein levels of E2F1 and SDMe. Actin was used as a loading control. (biological repeats: n = 4 for VCAN and METTL6, n = 3 for REV3L , and n = 8 for MDM1 ). See also . ( L ) A ChIP assay performed on WT E2F1 or E2F1 Cr HCT116 cells. Recruitment of E2F1 to the promoter regions of MDM1 and VCAN was tested. CDC6 acted as a positive control. Displayed is the mean percentage enrichment of input, with SD. Significance was calculated by Student’s t -test between the indicated sample pairs. An immunoblot is included to display input protein levels of E2F1 and SDMe. Actin was used as a loading control (biological repeats: n = 3 for MDM1, VCAN , and CDC6 ).
Article Snippet: Cleared extracts were then divided equally into tubes and mixed with 3 μg of E2F1 specific IgG (KH-95, sc-251, Santa Cruz, Dallas, USA; or clone G10, Argonaut Therapeutics, Oxford, UK) overnight.
Techniques: Generated, Western Blot, Control, Derivative Assay, Reverse Transcription Polymerase Chain Reaction, Positive Control